Thermodynamics of glasses: a first principle computation 
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Abstract 



We propose a first principle computation of the thermodynamics of simple 
fragile glasses starting from the two body interatomic potential. A replica for- 
mulation translates this problem into that of a gas of interacting molecules, 
each molecule being built of m atoms, and having a gyration radius (related 
to the cage size) which vanishes at zero temperature. We use a small cage 
expansion, valid at low temperatures, which allows to compute the cage size, 
the specific heat (which follows the Dulong and Petit law), and the configu- 
rational entropy. 
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Take a three dimensional classical system consisting of N particles, interacting by pairs 
through a short range potential. Very often this system will undergo, upon cooling or 
upon density increasing, a solidification into an amorphous solid state- the glass state. The 
conditions required for observing this glass phase is the avoidance of crystallisation, which 
can always be obtained through a fast enough quench (the meaning of 'fast' depends very 
much of the type of system) Jl|. There also exist potentials which naturally present some 
kind of frustration with respect to the crystalline structures and therefore solidify into glass 
states, even when cooled slowly- such is the case for instance of binary mixtures of hard 
spheres, soft spheres, or Lennard- Jones particles with appropriately different radii. These 
have been studied a lot in recent numerical simulations 0-§|. 

Our aim is to compute the thermodynamic properties of this glass phase, using the statis- 
tical mechanical approach, namely starting from the microscopic Hamiltonian. The general 
framework of our approach finds its roots in old ideas of Kauzman ||, Adam and Gibbs 0, 
which received a boost when Kirkpatrick, Thirumalai and Wolynes underlined the analogy 
between structural glasses and some generalized spin glasses ||10|| . In this framework, which 



should provide a good description of fragile glass-formers, the glass transition, measured 
from dynamical effects, is associated with an underlying thermodynamic transition at the 
Kauzman or Vogel-Fulcher temperature Tk- This ideal glass transition is the one which 
should be observed on infinitely long time scales 0. This transition is of an unusual type, 
since it presents two apparently contradictory features: 1) The order parameter is discon- 
tinuous at the transition: defining the order parameter as the inverse radius of the cage seen 
by each particle, it jumps discontinuously from in the liquid phase to a finite value in the 
glass phase. 2) The transition is continuous (second order) from the thermodynamical point 
of view: the free energy is continuous, and the transition is signalled by a discontinuity of 
the specific heat which jumps from its liquid value above T s to a value very close to that 
of a crystal phase below. These properties are indeed observed in generalized spin glasses 
The problem of the existence or not of a diverging correlation length is still an open 

one 
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This analogy is suggestive, but it also hides some very basic differences, like the fact that 
spin glasses have quenched disorder while structural glasses do not. The recent discovery of 



some generalized spin glass systems without quenched disorder |13[ has given credit to the 
idea that this analogy is not fortuitous. The problem was to turn this general idea into a 
consistent computational scheme allowing for some quantitative predictions. Important steps 
in this direction were invented in [0,0 , which showed the necessity of using several copies 
of the same system in order to define properly the glass phase. In a previous preliminary 
study, we used some of these ideas to estimate the glass temperature, arriving from the 
liquid phase p3fl . Here we concentrate instead on the properties of the glass phase itself, 
and particularly its properties at low temperatures. 
The Hamiltonian of our problem is simply given by: 

H= J2 v(x i -x j ) (1) 

l<i<j<N 

where the particles move in a volume V of a d-dimensional space, and v is an arbitrary 
short range potential. We shall take the thermodynamic limit N, V — > oo at fixed density 
p = N/V. For simplicity, we do not consider here the description of mixtures, which is 
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presumably an easy generalisation. The main obstacle to a study of the glass phase is the 
very description of the amorphous solid state. In principle one should give the average 
position of each atom in the solid, which requires an infinite amount of information. Had 
we known this information, we could have added to the Hamiltonian an infinitesimal but 
extensive pinning field which attracts each particle to its equilibrium position, sending N to 
infinity first, before taking the limit of zero pinning field. This is the usual way of identifying 
the phase transition. In order to get around the problem of the description of the amorphous 
solid phase, a simple method has been developed in the spin glass context- although one does 
not know the conjugate field, the system itself will know it, and the idea is to consider two 
copies (sometimes called 'replicas') of the system, with an infinitesimal extensive attraction. 
In the spin glass case this is a very nice method which allows to identify the transition 
temperature from the fact that the two replicas remain close to each other in the limit of 



vanishing coupling 24|25 



However this method is too naive and needs to be modified for the case of glasses. The 
reason has to do with the degeneracy of glass states. This property can be studied in detail 



in generalized spin glass mean field models fi5| , [14|| . For structural glasses, this is a conjecture 



which we shall make, on the basis of its agreement with the phenomenology of glasses ||. 
Let us assume that we can introduce a free energy functional F(p) which depends on the 
density p(x) and on the temperature. We suppose that at sufficiently low temperature this 
functional has many minima (i.e. the number of minima goes to infinity with the number 
iV of particles). Exactly at zero temperature these minima coincide with the mimima of 
the potential energy as function of the coordinates of the particles. Let us label them by 
an index a. To each of them we can associate a free energy F a and a free energy density 
fa — F a /N . The number of free energy minima with free energy density / is supposed to 
be exponentially large: 

Af(f,T,N)^exp(NE(f,T)), (2) 

where the function £ is called the complexity or the configurational entropy (it is the contri- 
bution to the entropy coming from the existence of an exponentially large number of locally 
stable configurations), which is not defined in the regions / > f max (T) or / < f m in(T), 
where J\f(f,T,N) = 0, and is supposed to go to zero at f m i n (T), as found in all existing 
models so far. In the low temperature region the total free energy of the system (fg) can 
be well approximated by the sum of the contributions to the free energy of each particular 
minimum: 

Z = e~ pNfs = e ~ mU - f maX d f e- N W-W> T V , (3) 

which shows that the minima which dominate the sum are those with a free energy density 
/* which minimizes the quantity $(/) — f — £(/, T)//3. The Kauzman temperature Tk 
is that below which the saddle point sticks at the minimum: /* = f m in(T). It is the 
only temperature at which there exists a thermodynamic singularity. Another characteristic 
temperature is the so called dynamical temperature T^: for Tp > T > Tk the free energy 
is still given the fluid solution with constant p and at the same time the free energy is 



also given by the sum over the non trivial minima [16,17|, and /* lies inside the interval 
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fmini fmax- The system may stay in one of the many possible minima. The entropy of the 
system is thus the sum of the entropy of a typical minimum and of £(/*, T), which is the 
contribution to the entropy coming from the exponentially large number of microscopical 
configurations. It is not obvious why this is equal to the liquid free energy. In this regime, 
the time to jump from one minimum to an other minimum is quite large: it is an activated 
process which is controlled by the height of the barriers which separate the different minima. 
The correlation time will become very large below To and for this reason To is called the 



dynamical transition point. It is also the mode-coupling transition temperature p0|-p2 
The real divergence of the correlation time appears at Tk- 

In order to get quantitative information on the behaviour of the system it is useful 
to consider the thermodynamics of m replicas which are constrained to stay in the same 



minimum |T7|]; this can be done introducing an extensive attraction among replicas which 
eventually goes to zero. In the same notation as before partition function is: 

J m 

which obviously coincides with the previous one for m = 1 In the limit where m — > 1 the 
corresponding partition function Z m is dominated by the correct saddle point /*, when the 
temperature is in the range Tk < T < To- For T < Tk, the saddle point /* sticks at 
/* — fmin(T) and the replicated free energy F m = — log(Z)/(/3m) is maximum at a value 
of m — m* smaller than one. One can use expression valid in the liquid phase (i.e. high 
temperature formulae) to evaluate the free energy F m at m < m*. We shall write down more 
explicit formulas in our case below. Notice that the 'replicas' which we introduce here play a 
slightly different role compared to the ones used in disordered systems: there is no quenched 
disorder here, and no need to average a logarithm of the partition function. 'Replicas' are 
introduced to handle the problem of the absence of description of the amorphous state. 
We do not know of any other procedure to characterize an amorphous solid state in the 
framework of equilibrium statistical mechanics. There is no 'zero replica' limit, but there 
is, as in disordered systems, an analytic continuation in the number of replicas. We shall 
see that this continuation looks rather innocuous. An alternative method is to introduce a 
real coupling of the system to another system which is thermalized [I6|; this has been used 
recently in order to study the glass phase PJISfl 

Let us turn to a more explicit implementation of these ideas. The original partition 
function, for N undistinguishable particles, is: 

Z =W\I ft(^) ex P ( ~P E - x i) ] ( 5 ) 

i=l \ l<i<j<N J 

We introduce m replicas of each particle, and compute Z m , in presence of an infinitesi- 
mal pinning field which is an attractive potential between them. This attractive potential 
(f)(xi, ...,xjsr) should not break the undistinguishability of all iV particles with the same 
replica index. We have found it convenient to use the attractive potential: 

exp(-(3(f>(xi,...,x N )) = — ex P -^EE^Kw -a 4(i))) ( 6 ) 

{TTaGSjv} \ i a,b J 
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where w is a small attractive potential which is short range (the range should be less than 
the typical interparticle distance in the solid phase), but its precise form is irrelevant. We 
then get a 'replicated' partition function: 



i „ N m I m N m \ 

• ~la=l V l<i<i<JVa=l i=l 8,6=1 ' 



A finite w gives rise to the formation of molecular bound states of m atoms. The ap- 
pearance of the glass states (T < Tp) is signaled by the fact that these molecules still exist 
in the limit lirn m _>i lirn TO _»o lini/v-+oo (notice the order of limits). According to the above 
discussion, the ideal glass transition (Tk) is detected from the existence of a maximum of 
the replicated free energy F m = — log(Z) / (/3m) at a value of m less than one. This is a well 
defined mathematical problem, which fully describes our general strategy for computing the 
thermodynamics of the glass state. Of course this cannot be done without resorting to some 
approximation schemes. We shall now develop one of them, a kind of harmonic expansion 
in the solid phase, but several other approximation schemes can be developed p6 |. 



We are interested in the regime of low temperatures where the molecules will have a 
small radius, justifying a quadratic expansion of v (we work here with a regular potential 
v(r), excluding hard cores). We thus write the partition function in terms of the center of 
mass and internal variables Zi,uf, with xf = + and J2a u i — 0, expand the energy to 
second order in u, and integrate over these quadratic fluctuations, leading to: 

m Nd/2/^ Nd ^~ l ) n ( m — 1 \ 
Z m = ^ JHdZiexpl-PmYtVizi-Zj) ~Tr log {f3M)\ (8) 

where the matrix M, of dimension Nd x Nd, is given by: 

= $ij E V nA Z i ~ Z k) ~ V^Zi - Zj) (9) 
k 

and v^ v {r) = d 2 v / dr^dr v (the indices \i and is, running from 1 to d, denote space directions). 
We have thus found an effective Hamiltonian for the centers of masses Zi of the molecules, 
which basically looks like the original problem at the effective temperature T* = l/(/3m), 
complicated by the contribution of vibration modes. We shall proceed by using a 'quenched 
approximation', i.e. neglecting the feedback of vibration modes onto the centers of masses. 
This approximation becomes exact close to the Kauzman temperature where m — > 1. The 
free energy is then: 

^ = -^ l ° g{m) ~ ^l^ l0g(27r) " l0gZ(n + ^(TrlogiPM)) (10) 

where the partition function Z(T*) is simply that of the usual monatomic liquid at the 
effective temperature T*, and the expectation value (.) is the Boltzmann expectation value 
at this same temperature. 

Let us notice that the condition for identifying the Kauzman temperature, 9 g^ n \ m =i — 0, 
reads in our harmonic approximation: 
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S liq = ^log(2vre) - \{Tr log {(5M)) (11) 

Su q is the entropy of the liquid at the effective temperature T e ff which equals T for m = 1. 
The right hand side of this equation is nothing but the entropy S so i of an harmonic solid 
with a matrix of second derivatives given by M. Thus we have found: 



d(3F n 



dm 



Suq — Ssol (12) 



m=l 



If Suq < S so i, one lies in the glass phase (T < T K ), while in the other case Su q > S so i, the 
temperature is greater than T K (and of course less than T D if the spectrum of M is positive). 
The complexity is then S c = Su q — S so i, as expected on general grounds [[17] . 



The harmonic expansion makes sense only if M has no negative eigenvalues, which is 
natural since it is intimately related to the vibration modes of the glass. Notice that here 
we cannot describe activated processes, and therefore we cannot see the tail of negative 
eigenvalues (with number decreasing as exp(— 1/T) at low temperatures), which is always 
present. It is known however that the fraction of negative eigenvalues of M becomes negligi- 
ble below the dynamical transition temperature Tjj [27]. So our harmonic expansion makes 
sense if the effective temperature T* is less than Tjj. 

Computing the spectrum of M is an interesting problem of random matrix theory, in 
a subtle case where the matrix elements are correlated. Some efforts have been devoted 
to such computations in the liquid phase where the eigenmodes are called instantaneous 



normal modes p7j. It might be possible to extend these approaches to our case. Here 
we shall rather propose a simple resummation scheme which should be reasonable at high 
densities-low temperatures. Considering first the diagonal elements of M, we notice that 
in this high density regime there are many neighbours to each point, and thus a good 
approximation is to neglect the fluctuations of these diagonal terms and substitute them by 
their average value. We thus write: 

(Tr log (J3M)) = Ndlog((3r ) + (Tr log - -v,„(z t - z k ))) (13) 

where 

1 f 

r ° = d J d dr 9{ r )^v(r) (14) 

and g(r) is the pair correlation in the liquid at the effective temperature T*. In principle 
the spectrum at this stage still depends on all the correlation functions of the liquid at T*, 
as can be seen from an expansion of (|13|) in powers of l/r . A simple 'chain' approximation 
involving only the pair correlation, consists of approximating in each term of order larger 
than 2 in this expansion the full correlation by a product of pair correlations: 



J dx\...dxp (^(xi, — , Xp) f^t 1 ^ 2 (^'i ^2) •••^Up—x/ip 1 •£p)'VfAptii ip^p "^1 

d d k 

W) 




r d d h / 1 \ p 

^ d -^Sw)-M k) -d m ) (15) 
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where the functions a and b are defined by: 



d d r g(r)v, v (ry kr = 5, u a{k) + - ^ b(k) (16) 

This chain approximation selects those contributions which survive in the high density limit, 
systematic corrections could probably be computed in the framework of the approach of 
|28fl , we leave this for future work. Here and in what follows, we have not written explicitly 



the density: we choose to work with density unity and vary the temperature (density and 
temperature variations are directly related in soft sphere systems onto which we focus below). 
The free energy within the chain approximation is: 

(m-1) r d Vflu (r) 2 



Am 



/^(r)E^ (17) 

J ,,,, ' n 



flu 

where the function L 3 is defined as: 

L 3 (x) = log(l - x) + x + x 2 /2 (18) 

We can thus compute the replicated free energy F m soley from the knowledge of the free 
energy and the pair correlation of the liquid at the effective temperature T* . We have done 
this computation in the case of soft spheres in three dimensions with v(r) = 1/r 12 , using 
the free energy and pair correlation function of the liquid given by the HNC approximation 
(obviously one could try to use better schemes of approximation for the liquid, depending 
on the form of v(r), in order to improve the results; our point here is not to try to get 
the most precise results, but to show the feasibility of a quantitative computation of glass 
properties using the simplest approximations). We find (always at density unity) a Kauzman 
temperature, obtained from the vanishing of (|T2D, which is ~ .194. When converted into 
the usual dimensionless parameter T = pT~ l l\ this gives r# — 1-51 which is close to the 
glass temperature (r ~ 1.6) observed in simulations with very fast cooling to avoid 
crystalization. Simulations done on binary mixtures (which do no crystallize) give a similar 
value for V . At the level of precision we have now reached, one will need both to do the 
theoretical computation for mixtures and also to perform careful simulations in order to 
disentangle the values of T K and T D . 
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FIGURES 




0.20 



FIG. 1. The full curve gives the inverse effective temperature 1/T* = (3m of the reference liquid 
as a function of the temperature T in the glass phase. The Kauzman temperature is Tx — .194. 
The dashed curves gives the square cage radius A divided by the temperature and multiplied by 
100. 

In fig. p] we show the values of the inverse effective temperature (1/T* = (3 * m) and of 
the square cage radius A, defined as A = §((xf) — (xi) 2 ). This square cage radius has been 
obtained by using in (^) as attractive potential: w{r) = r 2 /(4/3S), and differentiating the 
free energy: 



.4 



2 d{f3F) 
d{m-l)N 8(1/ S) 



(S = oo) 



(19) 



Notice that the effective temperature varies very little in the whole glass phase and 
remains close to the Kauzman temperature, while the square cage radius is nearly linear in 
temperature in the whole glass phase, which is natural since non harmonic effects have been 
neglected. The value of A at the Kauzman temperature is A ~ 2.5 1(F 3 . This corresponds 
to a typical lateral displacement of the particle in each direction of order \J~A ~ .05, which 
is .045 of the mean interparticle distance, a value which gives the correct order of magnitude 
for the Lindeman ratio. 

In fig. |2] we give the result for the specific heat in the glass phase versus the temperature. 
We see that it closely follows the Dulong-Petit law. This is the result that one should obtain 
since we study a solid phase in the classical framework. Notice that it is not at all trivial to 
derive this law from first principles in the glass phase. It is interesting to see it coming out 
naturally from our computations: although we are basically using the properties of the liquid 
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at the effective temperature T*, the fact that the optimal number of replicas m vanishes 
linearly with T at low temperatures naturally gives the Dulong-Petit law. 
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FIG. 2. The internal energy (full line) and the specific heat (dashed line) of the glass versus 
the temperature derived from our analytic computation. The specific heat value is close to the 
Dulong-Petit value of 3/2 and reaches this value in the zero temperature limit. 



From the knowledge of F m as a function of m, we can compute the configurational entropy 
as function of the free energy. In fig. |3|we plot the result for £(/) versus / at three different 
temperatures. We see that the curves are roughly parallel to each other, the main effect of 
the temperature changes being a shift in the / axis. 



9 



2.0 2.5 3.0 3.5 4.0 

f 

FIG. 3. The configurational entropy £(/) versus the free energy, computed analytically at 
temperatures T = .05, .1, .15, from right to left. 

As discussed above (see (0)), the value of the configurational entropy at equilibrium is 
zero for T < T#-. It becomes non zero above Tr-, where the saddle point in m is at m — 1. 
In fig. f| we plot the equilibrium configurational entropy versus the temperature. It will be 
interesting to try to compare it with numerical simulations. 




10 



FIG. 4. The equilibrium configurational entropy £ versus the temperature. It departs from 
zero above the Kauzman temperature. 

To summarize, we have developed a well defined scheme for the analytic study of the 
thermodynamics of the glass phase. The basic knowledge one needs is the detailed properties 
of the liquid (particularly the instantaneous normal modes) close to the glass transition. We 
have shown that an implementation of this scheme with rather simple approximations leads 
to very reasonable results. We hope to be able to refine these approximations in a near 
future in order to get very precise predictions. The extension of this approach to dynamical 
properties is also a fascinating perspective. 
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